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ABSTRACT 



Context. Massive stars play an important role in shaping the structure of galaxies. Infrared dark clouds (IRDCs), with their low 
temperatures and high densities, have been identified as the potential birthplaces of massive stars. In order to understand the formation 
processes of massive stars the physical and chemical conditions in infrared dark clouds have to be characterized. 
Aims. The goal of this paper is to investigate the chemical composition of a sample of southern infrared dark clouds. One important 
aspect of the observations is to check, if the molecular abuncances in IRDCs are similar to the low-mass pre-stellar cores, or whether 
they show signatures of more evolved evolutionary stages. 

Methods. We performed observations toward 15 IRDCs in the frequency range between 86 and 93 GHz using the 22-m Mopra radio 
telescope. In total, 13 molecular species comprising N2H+, '""CS, CH3CN, HC3N, HNC, HCO+, HCN, HNCO, CjH, SiO, H"CO+, 
H'^CN and CH3C2H were observed for all targets. Hence, we included in general species appropriate for elevated densities, where 
some of them trace the more quiescent gas while others are sensitive to more dynamical processes. 

Results. We detect HNC, HCO^ and HNC emission in all clouds and N2H^ in all IRDCs except one. In some clouds we detect 
SiO emission. Complicated shapes of the HCO^ emission line profile are found in all IRDCs. Both signatures indicates the presence 
of infall and outflow motions and beginning of star formation activity, at least in some parts of the IRDCs. Where possible, we 
calculate molecular abundances and make a comparison with previously obtained values for low-mass pre-stellar cores and high-mass 
protostellar objects (HMPOs). We show a tendency for IRDCs to have molecular abundances similar to low-mass pre-stellar cores 
rather than to HMPOs abundances on the scale of our single-dish observations. 

Key words. ISM: clouds, ISM: molecules. Radio lines: ISM, Stars: Formation 



1. Introduction 

Massive stars play an important role in determining physical, 
chemical, and morphological properties of galaxies. The last 
decade has seen considerable progress in th e und erstanding of 
massive star formation JZkinecker & Yor k3l2007l) . Since their 
detection bv iPerault et alj ( 11996 ) and E gan et alj ( Il998l) . cold 
and dense infrared dark clouds (IRDCs) appear to be the ideal 
sites to investigate the initial conditions for the process of mas- 
sive star formation. IRDCs have typical sizes between 1 and 
10 pc, masses from several hundreds to several thousands solar 
masses and H2 column densitie s between 2 and 10xlO^^_cm_^ 
(e.g. Rathborne et al. 2006; Vas vunina et al J 12009; Rag an et al] 
[2009 ). 

Apart from continuum observations, molecular line data 
have been used to characterize the proper ties of IRDCs. Th e 
first molecular line data were obtained by ICarev et all (Il998h . 
who detected H2CO in 10 clouds, thus confirming the pres- 
ence of dense gas. Using LVG modeling they estimated H2CO 
abundances of ~ 10"'*'. That is a factor of 50 lower in com- 
parison with low-density clouds and can be explained by accre- 
tion of gas-phase metals onto dust grains in the cold and dense 
IRDCs. Ammonia observations by Sridharan et al.. (,2005i) and 
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Imai et al.l (I2OO6I) allowed temperature determination of IRDCs 
in a range from 10 to 20 K. 

While earlier studies dSimon et al.ll2006l) suggested the asso- 
ciation of most IRDCs with the so-cal led Galactic molecu lar ring 
(galactocentric distance of 5 kpc), the I Jackson et alj i2008) study 
gave new evidence that the IRDC distribution in the first and 
fourth galactic quadrant more closely follows a galactic spiral 
arm (the Scutum-Centaurus arm). Since in normal spiral galax- 
ies, OB stars seem to form primarily in spiral arms, the associ- 
ation of IRDCs with a Milky Way spiral arm supports the idea 
that IRDCs are related to high-mass star formation. 

ISakai et alJ (l2008l) observed N2H+(l-0), HC3N(5-4), 
CCS(43-32), NH3(1,1), (2,2), (3,3) and CH30H(7-6) lines 
toward the massive clumps associated with IRDCs, to study the 
chemical conditions in IRDCs. Analysing the CCS and N2H^ 
abundance ratio, they conclude that infrared dark clouds are 
chemicaly more evolved than low-mass pre-stellar cores. 

An estimation of t he chemical evolutio nary status of IRDCs 
was performed also bv lGibson et al.l(l2009h . Using C"*0, CS and 
N2H^ abundances, and a chemical evolution code, they showed 
that cores where all three lines are detected appear to be chem- 
ically young (10^'^ < t < lO''''' years). Cores where no N2H^ 
emission is detected are suspected to be especially young (t < 
10^ years). This suggests that these regions may not have yet 
formed massive protostars. 
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A molecular survey by ISakai et aP (l2010l) toward 20 mas- 
sive clumps, including mid-IR bright and mid-IR dark sources, 
showed outflow activities and the possible presence of protostars 
in some clouds. Their line width analysis allowed them to recon- 
struct the distribution of the molecular species at the different 
evolutionary stages of massive clumps. 

While many characteristics of infrared dark clouds were de- 
termined during the last ten years, some of their properties are 
still not well known. Among the open questions is their chemical 
state. What is the chemistry in IRDCs? Is it really different from 
the chemistry in low-mass dark clouds? 

To enlarge the sample of well characterized IRDCs, we 
started a program to study gas and d ust properties of southern 
infrared dark clouds dVasvunina et al.t 2009. Paper I hereafter). A 
set of 15 clouds has been selected in the pre-Spitzer era by visual 
examination of the extinction contrast of the MSX 8.3 /im im- 
ages. In the meantime, the Spitzer satellite has since succeeded 
MSX and provided a far higher spatial resolution and sensitiv- 
ity. GLIMPSE mid-infrared images of our regions were retrieved 
from the Spitzer Archive. Continuum 1.2 mm maps were ob- 
tained with the SIMBA bolometer array at the SEST telescope. 
A 1.2 mm and 8 fim study of these southern IRDCs showed that 
these objects are not just distant Taurus-like clouds, but a distinct 
type of clouds with a clear trend to higher H2 column densities. 
It was found that the true peak column densities, extracted from 
millimeter data, exceed 3 xlO"-' cm"- (or 1 g cm"^), which has 
been proposed as a threshol d for high-mass star-forming clouds 
(iKrumholz & McKeell2008l) . 

This paper is our next step toward the understanding of the 
nature of the IRDCs. Here we present our investigations of the 
chemical composion of southern clouds. We perform molecular 
line observations in the 3 mm band with the Mopra single-dish 
radio telescope. Combining molecular line data and H2 column 
densities from the previous study we estimate molecular abun- 
dances and compare them with results for low-mass pre-stellar 
cores. Analysis of molecular lines provides not only information 
about chemistry, but also about physical processes in molecular 
clouds. For instance, the presence of SiO emission and HCO^ 
extended wings are evidences for outflow activity in a cloud. 
Specific line shapes can indicate infall motion. Detection or non- 
detection of some species help to determine the evolutionary sta- 
tus of our targets. 

The paper is organized in the following way. In Sect. 2 we 
describe our target selection, the selected molecular lines, obser- 
vational and technical details. In Sect. 3 we present the results of 
the qualitative analysis of obtained molecular line spectra, line 
parameters and abundance estimates. We discuss the obtained 
results in Sect. 4 and conclude in Sect. 5. 



2. Target and line selection, observations and data 
reduction 

2.1. Target selection 

A sample of 15 southern infrared dark clouds was selected by 
visual examination of the MSX 8.3 fim images for the pres- 
ence of high-contrast dark clouds. 1.2 mm continuum data 
were obtained for these objects with the SIMBA/SEST tele- 
scope. Together with millimeter data we used 8 /urn IRAC data 
from the Spitzer Galactic Legacy Infrared Mid -Plane Survey 
Extraordinaire (GLIMPSE, Benjami n et al.ll200 3^ to investigate 
the physical properties of the extinction and emission material 
(see Paper I). Based on our millimeter and mid-IR data and 
taking into account single-dish beam sizes at millimeter wave- 



lengths, we select several points for observations in every cloud. 
The criterion for selection was the presence of either a 1 .2 mm 
emission peak, a 8 jim emission peak, or a 8 fim extinction peak. 
Such a criterion let us cover IRDCs at the different evolutionary 
stages. Millimeter emission and mid-IR extinction indicate the 
presence of cold gas, typical for more quiescent regions. A typ- 
ical source of 8 jum emission usually contains a small infrared 
cluster and, therefore, is likely to be at a more advanced evo- 
lutionary stage than most dark regions in this IRAC band. The 
targets are listed in Table[T] 3-color Spitzer/GLIMPSE images of 
IRDCs together with the telescope and beam positions are pre- 
sented in the on-line material Fig. lA.lllA.3l 

2.2. Line selection 

In order to probe the dense and cold gas in IRDCs, we need 
appropriate tracers. The 3-mm band offers a large selection of 
molecular transitions. In particular, rotational transitions with 
low quantum numbers are accessible. In the following, we in- 
troduce the 3-mm lines we have chosen for our study. 

N2H^: This species is known to be a selective tracer of qui- 
escent gas (e.g.. ICaselli et al.ll2002h and is particularly suitable 
for studying the structure and kinematics of cold star-forming 
cores. The hyperfine structure allows to reliably measure the 
optical depth. N2H+ is known to be a "late depleter", thus it 
is not strongly affected by freeze out on grain surfaces (which 
prevents the use of the more common CO or CS transition). 
This makes it a robust tool for scrutinizing the highly structured 
interiors of massive star-forming regions where warm massive 
protostars might co-exist with younger cold and massive cloud 
cores. N2H^ has been - detected in both low- and high- mass 
pre-stellar cores and infrared dark clouds (e.g., Lee et alJl200ll : 
Pirogov et al."2003 UFuller et alJl2005l:lRagan et al.l l26oa 

'^CS: This molecule is a very good tracer of dense gas (e.g., 
iBronfman et al.ll 19961) due to its high di pole moment. It wa s de- 
tected in low-mass pre-stellar cores (e.g..l TafallaetalJ2002h and 
high-mass star forming regions (e.g., Jone s et alll2008h . In con- 
trast to N2H^, in cold dark clouds the CS emission vanishes to- 
ward the core center due to depletion, but can be used to trace 
layers surrounding the central cores. 

CH3CN: This mole cule is considered to be a tracer of warm 
and dense regions (e.g. jKalenskii et alJl2000l:lArava et alJl2005L 
I Vi til 12005 ). s Chemical models involving only gas phase reac- 
tions as well as models which take into account grain-surface 
chemistry, show that CH3CN is only detectable in an env iron- 
ment with elevated temperatures (see iPurcell et alj (|2006|) and 
references therein). It was detected in more evolved massive star- 
formation regions, like SgrB2 (Jones et al. 2008) orG305.2-i-0.2 
( i Walsh & Burton 2006.) , thus confirming the theoretical predic- 
tions. We include CH3CN in our "cold-gas-survey", because 
among infrared dark points we have quite a significant number 
of regions, where star-formation processes might have akeady 
started. 

HC3N: This molecule belongs to an important group of in- 
terstellar molecules - the cyanopolyynes, HC2„+i N. and is a valu- 
able tracer of physical conditions in molecular clouds. Since its 
first detection in space (Turner, 1971) HC3N has been found 
in every type of molecular cloud from giant molecular clouds 
associated with H II regions to circumstellar envelopes. It was 
shown that HC3N transitions have low optical depth and indi- 
cate the presence of denser gas than other high density tracers 
(Berginet al. 1996). 

HNC: This molecule is a commonly used tracer of dense gas 
in molecular clouds. The abundance ratio HCN/HNC strongly 



T. Vasyunina et al.: Chemistry in Infrared Dark Clouds. 



3 



depends on the temperature and in the case of the Orion molec- 
ular cloud decreases from 80 near the warm core to 5 on the 
colder edges (Goldsmith et al. 1986; Schilke et al. 1992|. The 
recent theoretical work by ISarrasin et al.i (12010 ) confirmed that 
this ratio should be around one for cold molecular clouds. This 
line is a triplet, but the spread of the hyperfine components is 
only 0.21 MHz or 0.4 km/s. 

HCO+ and H'^CO+: HCO+ is known to be a good tracer 
of th e dense gas especially of embedded molecular outflows 
(e.g. JCodella et alJl200UlHofner et al.ll2001h . HCO+ is an abun- 
dant molecule, with abundances especially enhanced around re- 
gions of higher fractional ionization. It is also enhanced by the 
presenc e of outflows w here shock-generated radiation fields are 
present dRawhngs et al. 2000, 2004). 

The emission from the H'^CO^ isotopologue is mostly op- 
tically thin and traces similar gas densities as HCO^. A com- 
parison between the generally optically thick HCO^ line and the 
optically thin H'^CO^ line yields information on the bulk mo- 
tion of the gas in the region. 

HCN and H'^CN: This molecule has been suggested as 
ubiquitous high density gas tracer. Moreover, the HCN molecule 
is known to be a good tracer of infall motions in low-mass star- 
forming regions. However, for the high-mass cores this can be 
different. Here HCN may become an unreliable infall trace r be- 
cause of a higher l evel of turbulence (iRedman et al.ll2008b and 
outflow signatures (IZhang et alJl20 07l). 

HNCO: IZinchenko etal] dJoOO) showed that HNCO inte- 
grated intensities correlate well with those of thermal SiO emis- 
sion in massive clouds. This can indicate a spatial coexistence of 
the two species and may hint to a common production mecha- 
nism, presumably based on shock chemistry. 

C2H: Chemical models predict that C2H is only well cen- 
tered on the sources when they are very young. At later stages it 
gets d estroyed in the centra l cores, but is replenished in the outer 
shells dBeuther et al.ll2008h . 

SiO: It can trace shocked gas potentially associated with en- 
ergetic young outflows. Hence, this line can reveal star formation 
activity e ven for cores w here no Spitzer/MIPS sources are appar- 
ent (e.g. jLo et alJl2007h . For a few of our IRDCs we have indi- 
cations from Spitzer/GLIMPSE imaging that shocked gas exists, 
often, but not exclusively, at the edges of the IRDCs. Tracers 
are the "green fuzzies" found in the 4.5 urn GLI MPSE channel 
(ICvganowskietal.ll2008t [Chambers et all I2009h and generally 
attributed to pure rotational IR lines of H2 (iDe Buizer & Vacca 
[2OIO). 

CH3C2H: A good tracer of de nse gas appropria te for 
early stages of star forming regions (iBergin et al.lll994 ). This 
molecule can be used as a good thermometer in a dense environ- 
ment. 

The covered molecular fines and transitions are summarized 
in Table |2] 

2.3. Mopra observations and data reduction 

Tfie observations were made witfi tfie 22-m Mopra radio tele- 
scope, operated by tfie Australia Telescope National Facility 
(ATNF), in tfie position switcfiing mode. In total we spent 7 
minutes on source and 7 minutes on tfie OFF position. Our tar- 
gets are dense molecular condensations witfiin larger molecular 
clouds witfi often widespread molecular emission. Tfierefore, we 
refrained from using one standard OFF position tfirow. Instead, 
OFF positions were cfiosen individually for every target region 
and were approximately 8-10' away from tfie source. 
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Fig. 1. Average spectrum. Obtained by averaging all spectra for 
every species witfi equal weigfit after sfiifting all of tfiem to tfie 
same reference velocity. Several weak detections are amplified 
by a factor of 10, and CH3CN is amplified by a factor of 20 for 
plotting. 



Tfie Mopra spectrometer (MOPS) off'ers zoom mode config- 
urations witfi tfie possibility to observe up to 16 sub-bands of 
138 MHz eacfi witfiin a total frequency range of 8.3 GHz. Tfiis 
set up delivers a velocity resolution of ~ 0. 1 1 km/s. 

Tfie observations were carried out on 9-11 May 2008 witfi 
tfie 3mm band receiver. We put tfie central frequency for tfie 8.3 
GHz block to 89270 MHz and tfius covered tfie range from 85 to 
93 GHz. In tfiis range we distributed 13 zoom windows, wfiicfi 
covered tfie lines listed in Table |2] 

System temperature measurements were performed every 30 
minutes and a pointing scan every fiour. Typical system temper- 
atures (measured witfi tfie common cfiopper-wfieel tecfinique) 
during tfie observations were 170-210 K. At Mopra observatory, 
SiO masers are used to correct tfie telescope pointing, giving a 
pointing accuracy better tfian 10". Tfie main beam of tfie tele- 
scope varies between 36" at 86 GHz and 33" at 115 GHz and 
tfie main beam e fficiency varies b etween 0.49 at 86 GHz and 
0.44 at 100 GHz dLadd et al.ll20"05h 

Mopra data are originally stored in RPFITS format. Using 
tfie ATNF Spectral line Analysis package (ASAP), we trans- 
formed tfiese raw data into ascii files wfiicfi were tfien fed into 
GILDAS for furtfier analysis. Tfie typical rms level in tfie ob- 
tained spectra is about 0.12 — 0.16 K. We give Icr errors in 
Tables ICTTIQ] 
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Table 1. List of observed IRDCs. 



Name R.A. Decl. distance" T'' N(H2)'' category^ 

(J2000.0) (J2000.0) (kpc) (K) *1022 cm-^ 



IRDC308.13-1 


13 37 01.582 


-62 44 34.01 


4.8 


13.5 


1.0 


Q 


IRDC308.13-2 


13 37 00.418 


-62 43 41.01 


4.0 


13.5 


1.0 


M 


IRDC308.13-3 


13 37 02.163 


-62 43 39.01 


4.0 


13.5 


1.0 


M 


IRDC309.13-1 


13 45 17.521 


-62 22 02.84 


3.9 


16.3 


0.6 


M 


IRDC309.13-2 


13 45 22.610 


-62 23 27.48 


3.9 


14.7 


0.5 


M 


IRDC309.13-3 


13 45 16.775 


-62 25 37.25 


3.9 


35.4 


0.4 


A 


IRDC309.37-1 


13 48 38.532 


-62 46 17.55 


3.4 


31.4 


0.9 


A 


IRDC309.37-2 


13 47 56.116 


-62 48 33.46 


3.4 


15.7 


0.8 


A 


IRDC309.37-3 


13 48 39.383 


-62 47 22.39 


3.4 


15.7 


0.7 


Q 


IRDC3 10.39-1 


13 56 01.359 


-62 14 18.29 


4.2 


27.4 


1.2 


A 


IRDC3 10.39-2 


13 56 00.759 


-62 13 59.80 


4.2 


27.4 


1.2 


A 


IRDC3 12.36-1 


14 11 27.752 


-61 29 27.18 


4.5 


13.5 


1.1 


A 


IRDC3 12.36-2 


14 11 56.773 


-61 29 25.78 


4.0 


14.4 


0.3 


Q 


IRDC3 13.72-1 


14 22 53.158 


-61 14 41.00 


3.3 


19.9 


0.4 


A 


IRDC3 13.72-2 


14 22 57.151 


-61 14 10.84 


3.3 


19.9 


0.4 


A 


IRDC313.72-3 


14 23 02.720 


-61 13 39.64 


3.3 


19.9 


0.3 


Q 


IRDC3 13.72-4 


14 23 04.533 


-61 14 46.00 


3.3 


19.9 


0.4 


Q 


IRDC316.72-1 


14 44 19.000 


-59 44 29.00 


2.7 


26.1 


1.2 


M 


IRDC316.76-1 


14 44 56.000 


-59 48 08.00 


2.7 


22.6 


4.1 


A 


IRDC3 16.72-2 


14 44 15.400 


-59 43 20.00 


2.7 


24.3 


1.3 


Q 


IRDC3 16.76-2 


14 45 00.500 


-59 48 44.00 


2.8 


23.2 


4.8 


A 


IRDC317.71-1 


14 51 06.905 


-59 16 11.03 


3.0 


15.6 


1.2 


Q 


IRDC317.71-2 


14 51 10.975 


-59 17 01.73 


3.0 


16.6 


3.5 


A 


IRDC317.71-3 


14 51 19.667 


-59 17 43.77 


3.2 


15.6 


0.6 


Q 


IRDC320.27-1 


15 07 56.251 


-57 54 32.11 


2.1 


15.3 


0.8 


Q 


IRDC320.27-2 


15 07 31.616 


-57 53 27.24 


2.2 


16.1 


0.5 


Q 
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IRDC321.73-1 


15 18 26.387 


-57 22 00.14 


2.2 


22.0 


1.0 


M 


IRDC321.73-2 


15 18 01.693 


-57 22 02.00 


2.2 


11.7 


1.7 


M 


IRDC32 1.73-3 


15 18 01.065 


-57 21 24.48 


2.1 


11.7 


1.7 


A 


IRDC013.90-1 


18 17 33.378 


-17 06 36.70 


2.5 


12.9 


2.6 


M 


IRDC013.90-2 


18 17 19.350 


-17 09 23.69 


2.4 


13.4 


1.1 


Q 


IRDC3 16.45-1 


1444 51.515 


-60 30 55.00 


3.1 


15.4 


0.7 


M 


IRDC3 16.45-2 


14 44 47.095 


-60 31 30.89 


3.1 


14.2 


0.5 


M 


IRDC318.15-1 


14 55 57.704 


-59 29 04.12 


3.0 


17.6 


0.5 


M 


IRDC318.15-2 


14 55 58.334 


-59 28 30.52 


2.9 


17.6 


0.4 


Q 


IRDC309.94-1 


13 50 54.970 


-61 44 21.00 


5.3* 


48.8 


5.2 


A 



Notes: Columns are name, right ascension, declination, distance, kinetic temperature, H2 column density and type according to mid-IR classifica- 
tion. 

" At first, to estimate the kinematic distances to our IRDCs, we used HCO^ line velocities (see Vasyunina et al. 2009). But more detailed investiga- 
tion showed that HCO^ lines, as a rule, have complex line shapes and the vlsr positions can be shifted up to 2 km/s in comparison with its optically 
thin isotopologue H'^CO^. We cannot use H'^CO^ for distance determination, since it is much weaker and we detected it not in all cases. Thus, 
for distance determination we decided to use N2H+, which is optically thin and distinguishable for all regions except IRDC309.37-2. Despite the 
significant difference in velocities (up to 2 km/s) between HCO^ and N2H^, the kinematic distances did not change drastically compared with the 

values in Paper I. 

* from lSaito et"^ ( 1200 ih 

Temperatures were derived based on (1,1) and (2,2) ammonia transitions observed with the the 64-m Parkes radio telescope (Linz et al. in prep.). 
To estimate H2 column densities in every point we used 1.2 mm data from SIMBA/SEST adopting the Mopra telescope beam size. 
"A" indicates "active" cores, "M" - "middle", "Q" - "quiescent". 



3. Results 

3.1. Qualitative assessment 

Figure[T]presents the average spectrum of an infrared dark cloud. 
To produce it we shifted the spectra of all our IRDCs to the same 
reference velocity and averaged spectra for each species with 
equal weight. The "average spectra" technique enables us to re- 
duce the noise level and recognize weak lines, which are not 
detectable in the single spectra. Spectra for every single point 
are available in the on-line material (see Fig. IB.HIB.19l l. For all 
IRDCs we have quite strong and clear detections of the HNC, 
HCO^ and HCN molecules. In all cases we see a non-gaussian 



shape of HCO^ with "shoulders" or double-peaked profiles. We 
will discuss the asymmetries of HCO^ in comparison with the 
optically thin H'^CO^ line in Sect. I3.3l in more detail. 

N2H^ is detected everywhere, except for IRDC309.37-2. As 
it was expected, we cannot distinguish all 7 hyperfine compo- 
nents because of the large line width. While one hyperfine com- 
ponent (I01-O12) constitutes a distinct line peak, the other six 
transitions merge into two satellites. 

HC3N, HNCO and C2H show rather weak emission and were 
detected in 18, 13 and 24 positions, respectively. 

According to th e previous studies (e.g. iRagan e t al. 200^ 
iJackson et al.ll2008h . CS is a common molecule in IRDCs. Still, 
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5 



Molecule 


Transition 


Rest frequency 
(GHz) 


A 

(*10"^s"') 


g« 


(K) 


Comments 


CH3C2H 


53-43 


85.442 600 


0.129778 


44 


1131 


Tracer of dense gas, 




52-42 


85.450 765 


0.170373 


22 


41.22 


good thermometer 




5i-4i 


85.455 665 


0.194760 


22 


19.53 






5o-4o 


85.457 299 


0.202908 


22 


12.30 




H'^CN 


li-O, 


86.338 735 








Tracer of dense gas, 




I2-O, 


86.340 167 


2.8 


9 


4.14 


infall motions 




lo-Oi 


86.342 256 










H"CO+ 


1-0 


86.754 330 


2.8 


3 


4.16 


Tracer of dense gas 


SiO 


2-1 


86.847 010 


2.0 


5 


6.25 


Trace shoked gas 


CiH 


1-0 3/2-1/2 F = 2-1 


87.316 925 


0.152757 


5 


4.19 


Tracer of early 
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we detected only very weak '^CS(2-1) emission in three clouds. 
The low detection rate of '^CS(2-1) in our case can be ex- 
plained by 50 times less abundant '^'CS in comparison with CS. 
Also depletion of CS (and '^CS) can play a significant role 
here (iBeuther & Henni ng 2009). Another molecule, CH3C2H, 
was detected only in IRDC316.76-1 and 1RDC317.71-2 and 
shows very weak emission in 1RDC316.72-1, 1RDC316.72-2 
and 1RDC3 16.76-2. 

Even in the most evolved regions 1RDC316.76-1 and 
IRDC3 16.76-2, where strong mid-lR and millimeter emission 
are present, we did not detect CH3CN as a typical hot core tracer 
Since CH3CN emission arises from a warm and compact region, 
one explanation of this non-detection can be the low spatial res- 
olution of the telescope and relatively low signal to noise ratio in 
our spectra. Using the "average spectra" technique we were able 
to reduce the noise level and to find very weak CH3CN emis- 
sion (see Fig.[TJ. This detection indicates that CH3CN is present 



at least in some parts of our clouds, but there is an abundance 
deficit for this molecule i n IRDCs. This result is in a greement 
with the previous work bv lBeuther & SridharanI (l2007h . 

Ac cording to ' Cvganowski et al.l (l2008h and [Chambers et al.l 
(l2009h it is more likely to see SiO, as an outflow tracer, in 
the clouds, where there is emission both at 4.5 and 24 //m. 
However, we detected SiO lines with extended line wings 
only in three regions with mid-lR emission: 1RDC313.72-1, 
1RDC313.72-2 and 1RDC316.76-1. Other potentially interest- 
ing regions with emission in mid-lR either show very weak SiO 
emission (1RDC317.71-2, IRDC309.94-1), or no SiO emission 
at an (1RDC309.37-1, IRDC3 10.39-2). The strongest SiO line 
was detected in the 1RDC32 1.73-1. Within the Mopra beam this 
region contains a very weak source at 24 yum and so weak emis- 
sion at GLIMPSE 4.5 ii m that it was not identifie d as extended 
green objects (EGOs) in ICvganowski et alj (I2OO8I) . We classify 
this object as "middle" (see below). Detection of SiO emission 
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Fig. 2. Line spectra for different mid-IR categories of IRDCs. (left) "Quiescent" cloud, (middle) "Middle" cloud, (right) "Active" 
cloud. 



in such source indicates that absence or extreme weakness of 

mid-IR emission does not mean that star formation processes 

I I II ^ 

are not taking place in a cloud (cf. Lo et al. 2007) and shows the 
necessity of spectral line observations to identify outflows and 
star formation activity in molecular clouds. 

Based on mid-infrared SPITZER data, [Chambers et al.l 

(|2009) subdivided infrared dark clouds cores in to "active" and 
"quiescent". Cores were classified as "active", if they showed 
emission both at 8 and 24 /im, and as "quiescent" if they con- 
tained neither IR emission signatures. We added a "middle" 
stage to this classification whenever detect emission at 24 mi- 
cron only and not at the shorter wavelength (see Table [TJ. In 
the Fig. [2] we presented the spectra f or typical clouds from ev- 
ery category. [Chambers et al.l (l2009h have shown that "active" 
cores have smaller sizes, higher densities and more pronounced 
water and methanol maser activity than the "quiescent" cores. 
However, from the qualitative analysis of the IVIopra spectra, we 
do not see clear molecular signatures of any of these three cate- 
gories of sources. 
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Fig. 3. Distribution of N^H^ line width for our I RDCs sample, 
more evolved region s from (Pirogov et al.ll2003h and low-mass 
pre-stellar cores from lLee et al. I d2001h . 



3.2. Line parameters 

Using the GILDAS/CLASS software, we fit the whole N2H+ 
structure, resulting in LSR velocities, full linewidths at half 
maximum, and optical depths (Tables ICTb . For the lines where 



there is no hyperfine splitting ('^CS, HC3N, HCO+, HNCO, SiO, 
H'^CO^), or it cannot be resolved (HNC), we estimate integrated 
and peak line intensities, LSR velocities and full linewidths at 
half maximum from gaussian fits (Tables lC2llC.8l ). For all lines 
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Fig. 4. Comparison between NaH^ line width and kinematic dis- 
tance to the objects for our IRDCs sample. 



Fig. 5. Number distribution of 5y - (V//co+-^/f"co+)/^^H"cc>+ 
for our sample (see Sect. 13.3b . Dashed lines indicate the iScr^v 
level. 



we estimated integrated areas, measured by summing the chan- 
nels between suitable velocity limits under the line (Table |C9] ). 
To compare our results with the line parameters of other IRDCs, 
low-mass starless cores and more evolved high-mass clouds we 
chose N2H^. This line has a relatively simple line shape - with- 
out extended shoulders and self-absorption features, and it is de- 
tected in 97% of the clouds. Moreover, N2H^ is widely detected 
in low- and high-mass pre-stellar and protostellar cores. 

We find that the line widths for our IRDCs vary in the range 
from 0.6 to 2.8 km/s and a re in agreem e nt wit h the results ob- 
tained for other IRDCs by Rag an et"an (l2006h . where the line 
width was 0.5-4.0 km/s with the mean value 2.2 km/s. In low- 
mass pre-stellar cores (Lee et al. 2001) strong NiH^ emission 
was detected in 75% of the objects. Lines in low-mass regions 
are narrow with line widths in the range of ~0. 2-0.4 km/s, and 
a mean value of 0.3 km/s, what corresponds to the thermal 
linewidth. High-mass proto stellar objects show N^H^ emiss ion 
in all selected targets (Pirog ov et al.ll2003l: [Fuller et al.ll2005l) . In 
comparison to low-mass cores, these targets have broader lines 
0.5-3.5 km/s. Figure [3]presents the distribution of the NaH^ line 
width for all three types of objects. 

Broader IRDC lines can indicate more turbulent processes 
within clouds. However, 36" beam, corresponds to 0.17 pc at 
1 kpc distance and 0.9 pc at 5 kpc. Thus, larger and larger vol- 
umes are covered by the same beam for regions further away and 
hence might offer a simple explanation for larger line widths. To 
check if the line widths become larger for more distant clouds, 
on the Fig. |4] we plot the line width against the kinematic dis- 
tance. For this comparison we use our N2H+ data from the 
Table IC.ll For our objects the correlation coefficient between 
the two parameters is only 0.2. Thus, large line widths of high- 
mass clouds cannot be explained exclusively by the large dis- 
tances. In the same vein, clumping, which can easily be hid- 
den within a 36" beam, can only partly explain such elevated 
line widths (iPillai et al. 2006). The linewidth - size r elationship 
also ca nnot explain high line width values in IRDCs. iHever et aP 
(1200 ih report that a clear line width - size relation is just ob- 
served in larger clouds complexes of > 10 pc in size. The IRDC 
clumps targeted by us are much smaller with typical sizes of 
around 1 pc. Furthermore, previous studies have akeady indi- 
cated that the line widths in high-mass star-forming clumps are 
much higher than would be predicted by the usual relations be- 



twee n line width and siz e determined from regions of lower mass 
(e.g. lShirleve"taDl200l . 

Comparing line parameters for "quiescent", "middle" and 
"active" regions, we can see that there is a trend to have larger 
line widths and higher integrated intensities in more evolved ob- 
jects. The mean line widths are 1.4, 1.7 and 2.2 km/s and the in- 
tegrated intensities are 4.7, 7.1 and 20.6 K km/s for "quiescent", 
"middle" and "active" sources respectively. These results are in 
good agreement with Chambers et al. (2009) who concluded that 
"active" sources are more evolved and present further evolution- 
ary stages in comparison with "quiescent" sources. 

We also utilized NaH^ to estimate kinematic distances to our 
objects (see Table [B- In case of IRDC309.37-2 the N2H+ line 
was not detected and we used HCO^ to estimate the distance. 
Since we observed several points within one cloud, distance de- 
termination helps us to confirm that clumps with associated mid- 
IR emission in IRDC309.37, IRDC310.39 and IRDC371.71 ai-e 
connected with the other "dark" parts of the respective clouds 
and are not just a projection effect. 

3.3. HCO^ and H^^C& line profiles 

As was already mentioned, the HCO^ lines show a non-gaussian 
profile in almost all clouds. In some cases even its isotopoloque 
H'^CO^ has a complicated line shape. Moreover, we see a sig- 
nificant shift, up to 2 km/s, between the peak positions of the 
optically thick HCO^ line and optically thin H'^CO^. A quan- 
titative estimate of this asymmetry could provide us with infor- 
mation about dynamical processes in the clouds. 

The most common way to extract line asymmetries is based 
on the comparison of optically thin a nd optically thick Un e 
positions: 6V = (VMck - V,hi„)/AVthin dMardones e t al.l ll997h . 
This estimate was wid ely used both for low- and high mass 
star formation regions. (Fuller et al." 2005t iPirogov et al.1 120031: 
Purcell et al. 2006; Lee et al. 2001). Using this criterion, we can 
subdivide our IRDCs in to "blue shifted" clouds with 6V < Q 
and "red shifted" clouds with 6V > 0. The blue excess could be 
due to infall motions and red excess expanding motions or out- 
ffow. Among the sources, where both HCO^ and H'^^CO^ were 
detected, we identified 19 "blue shifted" and 12 "red shifted" 
clouds (Fig.|5] Table IC. St . In order to exclude the sources, where 
differences between the V;,,- values for different species are dom- 
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inated by measurement errors, we have t o determine the t hresh- 
old value. For the low-mass starless cores Mard ones et al ] (fT997l) 
determined this threshold as 5cr, where cr is the typical dV error. 
In our case H'^CO^, which we use as optically thin component, 
for some sources is quite weak, which leads to the 5cr value equal 
0.5. Most of our clouds have \6V\ < 0.5, only 4 objects show a 
sufficient blue shift, and 6 a red shift. From these 10 candidates 
we also have to exclude sources with very weak and/or noisy 
spectra. After this selection we determine one infall candidate: 
IRDC3 13.72-4, and three sources show evidence of expand- 
ing motions: IRDC317.71-1, IRDC317.71-2 and IRDC317.71- 
3. However, this method allows us only to identify infall or out- 
flow candidates, but does not provide any quantitative values 
(e.g. infall speed). 

Another met hod to estimate line asymmetries quantitatively 
was described in lMvers et al.l (Il996h . They used an analytic two 
layer radiative transfer model to obtain theoretical spectra, com- 
pare them with the wide range of observed double- peak profiles 
and ex tract a formula for infall speed (see Eq. 9 from lMvers et al.l 
This method allows us to estimate the infall velocity in 
a cloud by measuring blue and red peak parameters of the op- 
tically thick line. Among all our targets we selected 3 clouds 
with clear red shoulders in HCO^ and non-complicated shape of 
H'^CO^ (see Fig. |6]l. The estimated infall velocities are around 
2 km/s. This value is much higher than the sound speed in such 
an environment - 0.3 km/s, or a t ypical infall spee d for low-mass 
pre-stellar cores 0.05-0.09 km/s dLee et al.ll20()TI) . 

3.4. Derived quantities: column densities and abundances 

To estimate column densities we assume LTE conditions and op- 
tically thin emission and use the following equation: 

N - 8a- gi 1 1 y 

^ g,e.p{-E,lkT,,) I Tmbdv, 

where A is the rest wavelength of the transition, A is the Einstein 
coefficient, g„ is the upper state degeneracy, }v{Tg^) and iyiTbg) 
are the equivalent Rayleigh-Jeans excitation and background 
temperatures, 0„,f is the partition function, and E/ is the energy 
of the lower level (iLee et al ] l2009h . In our calculations for g„, A 
and E/ we used values from The Cologne Database for Molecular 
Spectroscopy (CDMS) (Muller et al. 2001, 2005) (Table©. For 
the excitation temperatures T^x in all cases, except HCO^ and 
H'-'CO^ (see below), we assume ammonia kinetic temperatures 
(Linz et al. in prep.) (Table [T]i. We calculate the partition func- 
tion Qro, for every source by interpolating data from the CDMS 
for the particular source temperature T^x. 

Integrated intensities J T^i,dv for every species, were mea- 
sured by summing the channels between suitable velocity limits 
under the corresponding molecular line (Table |C9] l. 

The optically thin assumption is suitable for most of 
our lines: '^CS, HC3N, HNCO, C2H, SiO, Hi''CO+, H'^CN, 
CH3C2H. In case of optically thick emission (e.g. HCN, HCO^, 
HNC) Eq. [T] gives only a lower limit. However, for N2H^ and 
HCO^ we can make a better column density estimation. In the 
case of N2H^ we use the advantage of the hyperfine structure 
and estimate optical depths from the hyperfine components fit- 
ting. Then column densities for N2H^ can be calculated by mul- 
tiplying Eq. [T]by a factor of t/(1 - e '"). For HCO^ the pres- 
ence of an optically thin isotopoloque allows us to estimate 
m ore realistic v a lues f or HCO^. Using an algorithm described 
in If^ircell et alJ (l2006h and assuming that HCO^ is optically 



thick, H^^CO^ is optically thin and their relative abundance ratio 
X = [HCO*]/[H^'^CO^] equals 50, we calculated HCO+ column 
densities for the sources, where both HCO^ and H'^CO^ lines 
were detected. We do not apply the same algorithm for the pair 
HCN-H'^'CN, since H'^^CN is detected only in four objects. In 
spite of non-detection of CH3CN in any single spectra, the "av- 
erage spectra" technique let us recognize the very weak CH3CN 
(5o-4()) and CH3CN (5i-4i) lines (see Fig[T]). We use this av- 
erage spectra to estimate the integrated intensity J T^bdv and 
derive the upper limit column density for this species. 

To transfer column densities to abundances we use the 
N(species)/N(H2) ratio, where N(H2) values were estimated 
from the 1.2 mm SIMBA/SEST data (Paper I) adopting the 
Mopra telescope beam size (see Table [1]). In the case of CH3CN 
we use the average value for the molecular hydrogen column 
density of 1.2x10^^ cm"^ and derived upper limit for CH3CN 
abundances of 4. 0.xlO"". Abundances for other species are pre- 
sented in Table [3] 

For our column density and abundance estimates, systematic 
errors play the dominant role. One aspect is the assumption that 
all lines are optically thin. In case of moderately optically thick 
lines, the true column densities might be higher by a factor of 
2-3 for our objects. We also utilize kinetic temperatures derived 
from ammonia observations as excitation temperatures. This can 
give another factor of 2-5 to the estimated values. Systematic er- 
rors are also present in the molecular hydrogen column densities, 
derived from the millimeter continuum measurements, where we 
made assumptions about the dust model and equal gas and dust 
temperatures. All these assumptions give rise to a systematic er- 
ror of an order of magnitude in molecular abundances which we 
obtained. 

3.5. Are tlie objects really cold? 

To estimate kinetic temperatures in our clouds we observed the 
ammonia (1,1) and (2,2) inversion transitions with the Parkes ra- 
dio telescope. These transitions are known as a good thermome- 
ter for cold, dense gas. We derived kinetic temperatures for most 
IRDCs between 10 and 30 K (Linz et al. in prep.). Ammonia 
observations show the presence of cold gas in IRDCs. However, 
the Parkes telescope beam is 72" or approximately 1 pc at the 
typical IRDCs distances. With the Mopra telescope (telescope 
beam 36") we can trace a smaller volume and use other molec- 
ular lines to identify warmer regions. 

One way is to consider the HCN/HNC abundance ratio. This 
ratio strongly depends on the temperature and in the case of the 
Orion molecular cloud decreases from 80 near the warm core 
to 5 on t he colder edge s (fGoldsmith et al. 1986; Schilke et aTj 
1992; Sa rrasin et al.l|20i 0). For our IRDC sample the estimated 
HCN/HNC abundance ratio values are not higher than 3. This 
confirms, that we deal with quite cold clouds. 

Beside ammonia, symmetric top molecules like CH3C2H can 
be used as probes of kinetic temperature in molecular clouds 
(iBerginet al.lll994 . We detect CH3C2H in 5 targets. In two ob- 
jects the 52-4-2, 5i-4i and 5o-4() components are clearly seen. 
This allows us to estimate rotational temp eratures in these re- 
gions, using the algorithm given in lBerginetal.l(ll994 . Derived 
temperatures equal 17 K for IRDC317.71-2 and 38 K for 
IRDC316.76-1. The temperature for IRDC317.71-2 is similar to 
the value we have from ammonia observations, whereas the tem- 
perature for IRDC316.76-1 is 15 K higher This indicates that 
in some clouds a warmer component can be present in addition 
to cold gas. In the case where there is also a warm gas com- 
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Fig. 6. Samples of the HCO^ and H'^CO^ line profiles. The vertical lines indicate the Wisr positions measured from H'^CO^ Une 
profiles. 
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Fig. 7. (left) Comparis on of the mole cular abund ances of the I RDCs with low-mass pre-stellar cores from 'Padov ani et aH (l2009l) . 
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ponent and the real gas temperature in some targets is higher, 
we underestimated molecular abundances. The abundances for 
IRDC316.76-1 with T=38 K become a factor of 1.2-1.7 higher. 
However, taking into account all our assumptions and the still 
large telescope beam, this difference is still within the error limit. 

3.6. Special objects 

In this subsection we present infrared dark clouds, where line 
spectra show some unique features. 

IRDC317.71 One of the remarkable features in this object 
is the ratios of HCN hyperfine components, which are differ- 
ent from other IRDCs (see Fig |2]). In contrast to N2H^, the 
HCN hyperfine com ponent intensity ratio is not a constant value. 
lAfonso et al.l (Il998l) show for low-mass dark clouds, that HCN 
hyperfine components can have different intensity ratios in dif- 
ferent sources. For most of our IRDCs we detected that the ratios 
F=2-l to F=0-1 and F=2-l to F=l-1 components are different 



but larger than 1. For IRDC317.71-1 these ratios are 0.8 and 0.4 
and for IRDC317.71-2 are 0.5 and 0.2 respectively. 

Another interesting feature is the large shift between 
H'^CO^ and HCO^ emission peaks. H'^CO^ emission peak is 
shifted to the blue side in comparison with HCO^ and corre- 
sponds to the extinction in the optically thick component. Also 
IRDC317.71-2 is one of the five locations where CH3C2H was 
detected. As we can see from the previous section, the CH3C2H 
excitation analysis shows the presence of a warmer gas compo- 
nent thus indicating a more evolved star forming region. 

Most probably the HCN spectrum can be explained by fore- 
ground absorption in low-excitation molecular gas. Similar hy- 
perfine anomalies are frequently seen in cold dark clouds where 
they are explained in this way (absorption in the envelope). The 
same effect is probably seen in HCO"^. The HCO^ line asym- 
metry apparently indicates the velocity shift of the foreground 
absorption. If attributed to the envelope, this shift indicates an 
expansion. In addition, HCO^ perhaps traces a contribution from 
an outflow. 
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IRDC317.71-2 corresponds to the bright emission source 
in all mid-IR Spitzer/GLIMPSE bands. This, together with all 
spectral features, which we discuss above, indicate ongoing star- 
formation processes in 1RDC317.71-2. Cross correlation of the 
mid-lR data and near-lR 2MASS K-band images show that the 
bright mid-lR source hides a young stellar cluster. 

IRDC321.73 In 1RDC32 1.73-1 (Fig|2]i, in spite of no emis- 
sion in GLIMPSE bands and only a very weak source at 24 fim, 
we detected clear SiO emission. This indicates the presence of 
a shock and, probably, outflow activity in this region. The pres- 
ence of the red HCO^ shoulder may be an evidence of infall 
motion. The other two points in this cloud (1RDC32 1.73-2 and 
IRDC32 1.73-3) show a more complicated than average shape of 
HCO^. The rest of the lines, even optically thin like N2H^, also 
have non-gaussian shapes. This can be explained by the presence 
of several kinematic components. 



Table 3. Abundances of the IRDCs. a(b) denotes a* 10* 
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" "A" indicates "active" cores, "M" - "middle", "Q" - "quiescent". 
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4. Discussion 



ICarev et al.l (Il998h and lPillai et~an (l2006h mention a difference 
in the ammonia and formaldehyde abundances and hence a pos- 
sible chemical difference between low-mass pre-stellar cores and 
IRDC s. To study the chemical conditions in IRDCs, Sakai et al. 
(l2008h observed N2H+(l-0), HC3N(5-4), CCS(43-32), NH3(1,1), 
(2,2), (3,3) and CH30H(7-6) lines toward the massive clumps 
associated with IRDCs. They estimated the CCS and N2H+ 
abundance ratio and used low-mass prestellar core abundances 
from Benson et al. ( 1998) for the comparison. The result of this 
analysis showed that IRDCs have lower N(CCS)/N(N2H+) ra- 
tio values and hence, they can be chemi cally more evolved 
than low-mass pre-stellar cores. However, ISakai et al] (l2008l) 
used particular low-mass objects for the comparison and did not 
take into account the pre-stellar cores L1512 and L63, where 
N(CCS)/N(N2H^) ratio is low and is in agreement with values 
obtained for their IRDC sample. 

Here we want to determine if the molecular abundances in 
IRDCs are similar to the low-mass pre-stellar cores, or whether 
they show signatures of more evolved evolutionary stages. To 
answer this question we compare our IRDC abundances with the 
data available in the literature for low-mass starless clouds and 
more evolved high-mass objects. We analyse the dispersion of 
the particular abundance values and the difference between mean 
abundances for all available species. As an additional criteria, we 
compare values of the HNC/HCN abundance ratio for different 
types of objects. 



4.1. Abundance comparison with low-mass starless cores 

To make a comparison wi th low-mass cold clouds, we used 
molec u lar abundances f:r om Padovani et al. (l2009l). Tafalla et aP 
(I2OO6I) . lJ0rpensen et al.1 (12004') and ' Ohishi et al.l (1 19921) (Fig. 171 
left panel). Comparing the dispersion of the abundance values 
we find that for some species, like HC3N, there is a large spread 
in abundance for the low-mass sample, but higher and lower val- 
ues occur with regard to the IRDC values. In case of HNCO 
we compute abundances for low-mass starless cores using new 
results by Marcelino et al. (2009) and adop t ing th e H2 column 
densities from Table 3 of IMarcelino et al.l (l2010h . That gives 
us quite low abundance values for Bl, L1527-b, and L1544. 



However, cores like L183 and TMC-1 reach HNCO abundances 
which are overlapping with IRDC values and just a factor of 2 - 
3 lower than the IRDC average abundance. The same situation, 
in principle, is found for HCO^, where an overlap between the 
low-mass core and IRDC HCO^ abundances exist for the ob- 
Jects LI 83 and TM C-1 (both attaining 8.0 x 10"^ according to 
IOhishietal.l([T992l) ). There is a tendency for higher abundances 
in IRDCs compared to low- mass cores for N2H^. However, also 
here we found with LI 544 a prominent example of a low-mass 
prestellar core with a rather high N2H^ abundance of 5.0 x 10"^ 
(J0rgensen et al. 2004). 

The difference between mean abundances of IRDCs and low- 
mass objects is a factor of 2-3 for all species and around 5 in 
the extreme case of HNCO (see Table O. Taking into account a 
large spread in abundances for every molecule within our IRDCs 
sample and systematic errors (see Sect. 13.4b . such differences in 
average abundances between the low-mass pre-stellar cores and 
IRDCs cannot be used to claim a clear and unequivocal chemical 
distinction between the two groups of objects. 



4.2. Abundance comparison with high-mass protostellar 
objects 

The next class of objects, which is interesting for the compari- 
son with our results, is the class of more evolved high-mass pro- 
tostellar objects (HMPOs). HMPOs are bright at mid- and far- 
IR wavelengths and characterized by higher gas temperatures, 
compared to IRDCs. It is assumed that these objects present 
one of the early stages of high-mass star formation, where the 
cent ral protosta r has a mass > 8 Mp and is still accreting, (e.g. 
Sri dharan et a n|2002; Beut her et al)|2 007). For the comparison 
we use data from Blake et al] (11987'). Helmich & van Dishoeckl 
(ll997h . |Pirogov et al.1 (l200'3r and Zinchenko et al. (2000, 2009^ 
HMPOs and IRDCs have different physical conditions (e.g. tem- 
peratures) and, hence, are expected to show different chemistry. 

From our comparison we can see, that high-mass proto- 
stellar objects show higher mean HC3N and HCN abundances, 
since some HMPOs have extremely high abundances of these 
species. C2H, on the contrary, has the mean abundance a fac- 
tor of 20 lower in the case of HMPOs. We did not detect 
CH3CN emission in any single IRDCs spectra. However, us- 
ing the "average spectra" technique we reduce the noise level 
and find very weak CH3CN emission and estimate its up- 
per abundance limit. This upper limit is lower in compari- 
son with the HMPOs a bundances from single-dish studies by 
dHelmich & van Dishoe ck 1997) and (Blake etal. 19 8^) (see 
Fig. Q right panel). Previous theoretical ( Nomu ra & Millaii 
2004) and observational (Beuther & Sridharan 2007) studies of 
high-mass star forming regions showed that a low amount of 
CH3CN is typical for the earliest evolutionary stages. 

The average abundances of HC3N, HCO+, HCN and C2H for 
IRDCs and HMPOs differ by a factor of 16-75. These differ- 
ences are significant, even when considering all of the included 
assumptions and uncertainties. 

4.3. HCN/HNC abundance ratio 

Another criterion that we use to determine the chemical status 
of IRDCs is the HCN/H NC abundanc e ratio. F or the low-mass 
preste llar c ores from 'Tafalla et"an (l2006h . Ij0rgensen et al.l 
(|2004 and lOhishi et al., (.1992) this ratio is < 1. Values of 
the HCN/HNC abundanc e ratio for HMPOs are higher and 
reach values up to 13 (^see iHeknich & van D ishoeck' ("1997^ and 
Blake et al. ( 1987)) a nd are around 80 in th e extreme regions 
Uke Orion (Gol dsmith et all [19861 ISchilke et al] Il992l) . We 
find the HCN/HNC abundance ratios for our IRDCs ~ 1. This 
value is in agree ment with the t heoretical predictions for the 
cold clouds ( Sarrasin et al. l l2OT0i) and values for the low-mass 
pre-stellar cores. 

Considering the comparison between IRDCs, low-mass star- 
less cores and HMPOs and assuming that the HCN/HNC abun- 
dance ratio strongly depen ds on temperature and is enhanced in 
active star-forming cores ( iGoldsmith et alj 119861: ISchiUce et aTl 
Il992h . our results support the idea that IRDCs present rather low 
mass star-less core chemistry than HMPOs chemistry. However, 
the number of low-mass pre-stellar cores and HMPOs, which 
we used for the analysis, may be not sufficient in a strict statis- 
tical sense. Therefore, to make more solid statement about the 
evolutionary status of IRDCs and be more confident, we need to 
extend the line samples, add sulfur bearing species, and perform 
further studies, including chemical modeling. 

Beside the comparison with low- and high-mass cloud abun- 
dances, we also compare molecular abundances of "quiescent". 
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"middle" and "active" regions within our IRDC sam ple. In ac- 
cordance with recent results bv Battersbv e t al.l (l20IOh . we do not 
detect any significant difference among these three categories. 

5. Conclusions 

In this paper, we present 3 mm molecular line observations with 
the 22-m Mopra radio telescope. In total 1 3 molecular lines were 
observed for all IRDCs. The results of our study can be summa- 
rized as follows: 

1 . Using H2 column densities from the previous investigation, 
we estimate molecular abundances of all species. We show 
that there is a tendency for the IRDCs to have molecular 
abundances similar to the low-mass pre-stellar core rather 
than to the HMPOs abundances. However, the derived abun- 
dances come with uncertainties of around one order of mag- 
nitude. Furthermore, also the comparison abundances for 
low-mass cores and HMPOs can be affected by consider- 
able uncertainties, especially, when abundances have been 
computed in the literature by combining heterogeneous data 
sets (regarding beam sizes etc.). To make more solid state- 
ments about the evolutionary status of IRDCs therefore calls 
for subsequent systematical studieS; 

2. According to the classification of Chambers et al ](|2009), we 
subdivided our clouds to "quiescent" and "active" and added 
"middle" class to them. We have found a trend for more 
evolved regions to have higher line widths and integrated in- 
tensities. However, we do not detect clear evidence of differ- 
ent chemistry in these three groups. 

3. Comparison of the line width and integrated intensities of the 
IRDCs and low-mass dark clouds show several times higher 
values for IRDCs. Broader and more intense lines mean that 
in IRDCs we have more turbulent conditions in comparison 
with low-mass clouds. 

4. We detect the SiO emission in some clouds and compli- 
cated shapes of the HCO^ emission line profile in all IRDCs, 
which indicates the presence of infall and outflow motions 
and the beginning of star formation activity, at least in some 
parts of the IRDCs. 

5. The analysis of the two available CH3C2H excitation di- 
agrams and detection of the very weak CH3CN (5o-4()) 
and CH3CN (5i-4i) lines on the "average spectra" indicate 
the presence of a warm gas component in some IRDCs. 
However, these warm regions are compact and cannot be re- 
solved with single-dish observations. 
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Appendix A: 3-color Spitzer/Glimpse images. 



T. Vasyunina et al.: Chemistry in Infrared Dark Clouds., Online Material p 2 

IRDC013.90 IRDC308.13 




RIGHT ASCENSION RIGHT ASCENSION 

IRDC309.13 IRDC309.37 




RIGHT ASCENSION RIGHT ASCENSION 

IRDC30994 IRDC310.39 



S 4+' 30" 




13 51 08 04' 



50 56 
RIGHT ASCENSION 




RIGHT ASCENSION 



Fig. A.l. 3-color Spitzer/GLIMPSE image of the Infrared Dark Cloud, where 3.6 jum is blue, 4.5//m is green and 8fim is red. Circles 
mark observed positions and show the beam size. 
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Fig. A.2. 3-color Spitzer/GLIMPSE image of the Infrared Dark Cloud, where 3.6 jum is blue, 4.5//m is green and 8/im is red. Circles 
mark observed positions and show the beam size. 
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Fig. A.3. 3-color Spitzer/GLIMPSE image of the Infrared Dark Cloud, where 3.6 //m is blue, 4.5jum is green and 8yum is red. Circles 
mark observed positions and show the beam size. 
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Appendix B: Spectra images. 
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Fig. B.l. Observed lines 
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Fig. B.2. Observed lines 
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Fig. B.3. Observed lines 
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Fig. B.4. Observed lines 
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Fig. B.5. Observed lines 
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Fig. B.6. Observed lines 
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Fig. B.7. Observed lines 
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Fig. B.8. Observed lines 
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Fig. B.9. Observed lines 
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Fig. B.IO. Observed lines 
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Fig. B.ll. Observed lines 
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Fig. B.12. Observed lines 
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Fig. B.13. Observed lines 
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Fig. B.14. Observed lines 
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Fig. B.15. Observed lines 
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Fig. B.16. Observed lines 
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Fig. B.17. Observed lines 
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Fig. B.18. Observed lines 
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Fig. B.19. Observed lines 
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Appendix C: Line parameters. 



T. Vasyunina et al.: Chemistry in Infrared Dark Clouds., Online Material p 26 
Table C.l. N2H'^ line parameters. 



Name 


V,,„- 


A V 




1 (T rms 




(kni/s) 


(km/s) 




(K) 






3.0 


(0.24) 




n 1 5 




-47 9 r)4'l 


0.7 


(0.11) 




0.16 


TRnr'^fiR 1 ^ ^ 




0.9 


(0.08) 




n 14. 






1.8 


(0.16) 




n 1 5 


LlS.LJ\^J\jy . ID A 




2.3 


(0.12) 




n 13 


ii^u\^j\jy , i-j J 


-4S CO 04^ 


0.9 


(0.13) 




n 13 




49 CO 09^ 


2.4 


(0.00) 


9 CO 0071 






-49 S CO OSI 


2.4 


(0.13) 




n 13 




-S9 5 CO on 


2.6 


(0.06) 


Q CO 9C^01 


0.12 




-S9 s cn 0^^ 


1.8 


(0.12) 


1 ft CO 7441 

1 .U [\J. I'-r'-t) 


n 13 




-SO 7 (0 OW 


2.4 


(0.08) 




0.12 




-SI WO 09^ 


0.7 


(0.04) 




n 1 1 


TRDr'^l 79-1 


-AA 1 (0 01\ 


2.3 


(0.09) 




0.12 


TRDr'^l 7?-? 


-AA CO n 

-'-T'-T.yj \ \J.\J L } 


2.0 


(0.07) 




0.12 




_A'\ S CO 09^ 


1.7 


(0.09) 




0.12 


TRDP'^I'^ 79-4 


-41 9 CO 0^^ 


2.1 


(0.10) 




0.12 


TRDP'^lfi 79 1 


CO 01 ^ 


2.9 


(0.04) 


1 Q CO 9 1 41 


n 19 


TRDP'^lfi 76-1 

IIVL/V^ J 1 \J, 1 \J 1 


_'XQ 7 CO 001 


2.8 


(0.00) 


8 CO 01 S1 


0.12 




s CO 09"^ 


2.2 


(0.06) 


9 ^ CO ^^791 

Z.J l^V./.J / z^ 


0.12 




-40 s CO on 


2.5 


(0.03) 


1 4 CO 1 C^CS1 

i .'-T 1 UU ) 


0.12 


TRDr"^ 17 71 1 




2.1 


(0.06) 


1 R CO ■^RQ1 


fl 1 9 


TRr)r'^17 71-9 


-A% (S CO 001 


2.4 


(0.00) 


Q CO 1 9Q1 
\).y 1 £.y ) 


0.12 


TRr)r'^17 71-^ 


-4^ q CO on 
-'-TD.y \\).\) i ^ 


1.8 


(0.14) 


1 4 CO 1 SCSI 


0.12 


TRr)r^90 97-1 


-^1 Q CO 041 


1.0 


(0.09) 


A'\ (\ Q1 81 


n 13 


TRnr'^90 97-9 


-'^9 4 CO 091 


0.5 


(0.06) 




n 13 




-X> CO OC^I 


0.5 


(0.06) 




n 07 




JV) CO oil 


1.7 


(0.06) 


1 4 CO ^Rfil 


n 13 


IRDC321.73-2 


-33.3 (0.03) 


1.5 


(0.13) 


4.5 (1.172) 


0.12 


IRDC321.73-3 


-31.8 (0.03) 


2.9 


(0.07) 




0.12 


IRDC013.90-1 


23.1 (0.03) 


1.6 


(0.10) 


1.8 (0.635) 


0.12 


IRDC013.90-2 


22.5 (0.03) 


1.3 


(0.10) 


2.6 (0.967) 


0.11 


IRDC3 16.45-1 


-43.5 (0.02) 


1.4 


(0.08) 


5.4 (0.994) 


0.13 


IRDC3 16.45-2 


-43.4 (0.03) 


1.6 


(0.10) 


4.3 (1.063) 


0.14 


IRDC318.15-1 


-42.6 (0.03) 


1.4 


(0.11) 


3.4 (1.273) 


0.12 


IRDC318.15-2 


-42.0 (0.02) 


1.0 


(0.08) 


2.7 (1.093) 


0.12 


IRDC309.94-1 


-59.2 (0.01) 


2.5 


(0.05) 


1.1 (0.246) 


0.16 



Notes: Columns are name, LSR velocity, full linewidth at half maximum, main group optical depth, Icr rms values. 



Table C.2. '^CS line parameters. 



Name 




V,,„- 


A V 


T 


1(T rms 




(K km/s) 


(km/s) 


(km/s) 


(K) 


(K) 


IRDC313.72-2 


0.3 (0.02) 


-44.1 (0.1) 


1.3 (0.2) 


0.2 


0.06 


IRDC3 16.76-1 


0.4 (0.01) 


-39.8(0.1) 


1.0 (0.0) 


0.4 


0.06 


IRDC3 16.76-2 


0.5 (0.02) 


-40.6 (0.1) 


1.9 (0.2) 


0.3 


0.06 



Notes: Columns are name, integrated area obtained from the gaussian fits, LSR velocity, full linewidth at half maximum, main beam brightness 
temperature, Icr rms values. 
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Table C.3. HC3N line parameters. 



Name 


iTmr dv 




A V 


T 


Icr rms 




(Kkm/s) 


(krn/s) 


(km/s) 


(K) 


(K) 


IRDC309.37-1 


2.1 (0.09) 


-42.6(0.10) 


2.3 (0.2) 


0.8 


0.42 


IRDC3 10.39-1 


1.0 (0.03) 


-52.8 (0.06) 


2.4 (0.1) 


0.4 


0.13 


1RDC3 10.39-2 


0.4 (0.02) 


-52.5 (0.07) 


1.5 (0.2) 


0.2 


0.13 


IRDC312.36-1 


0.3 (0.02) 


-50.7 (0.09) 


1.4 (0.2) 


0.2 


0.12 


IRDC313.72-1 


1.0 (0.03) 


-44.0 (0.09) 


3.0 (0.2) 


0.3 


0.13 


IRDC313.72-2 


1.0 (0.02) 


-44.2 (0.05) 


2.1 (0.1) 


0.4 


0.12 


IRDC3 16.72-1 


1.0 (0.02) 


-39.0 (0.07) 


2.6 (0.1) 


0.3 


0.12 


IRDC3 16.76-1 


4.2 (0.03) 


-39.9 (0.02) 


3.0 (0.1) 


1.2 


0.12 


1RDC3 16.72-2 


0.5 (0.02) 


-39.5 (0.06) 


1.8(0.1) 


0.3 


0.12 


1RDC3 16.76-2 


2.3 (0.02) 


-40.8 (0.02) 


2.2 (0.0) 


0.9 


0.12 


1RDC317.71-1 


0.5 (0.02) 


-44.2 (0.12) 


2.3 (0.2) 


0.2 


0.13 


TPrtP'^ 17 71 9 


J.\j yJ.V^} 


AA 1 /^n on 




1 A 


1 ^ 
U. iO 


IRDC320.27-1 


0.2 (0.01) 


-31.9 (0.05) 


0.8(0.1) 


0.3 


0.13 


IRDC321.73-1 


1.3 (0.02) 


-32.3 (0.03) 


1.8 (0.1) 


0.6 


0.13 


IRDC321.73-2 


1.0 (0.02) 


-33.7 (0.03) 


1.4(0.1) 


0.6 


0.13 


IRDC32 1.73-3 


0.7 (0.04) 


-31.3 (0.10) 


2.4 (0.5) 


0.2 


0.13 


IRDC013.90-1 


0.4 (0.02) 


22.8 (0.05) 


1.2(0.1) 


0.3 


0.12 


IRDC316.45-1 


0.6 (0.03) 


-43.4 (0.09) 


2.2 (0.2) 


0.2 


0.15 


IRDC309.94-1 


2.7 (0.04) 


-59.5 (0.04) 


2.9(0.1) 


0.8 


0.16 



Notes: Columns are name, integrated area obtained from the gaussian fits, LSR velocity, full linewidth at half maximum, main beam brightness 
temperature, Icr rms values. 



T. Vasyunina et al.: Chemistry in Infrared Dark Clouds., Online Material p 28 
Table C.4. HNC line parameters. 



Name 


J Tmb dv 
(Kkm/s) 


(km/s) 




AV 
(kni/s) 


T mb 

(K) 


Icr rms 
(K) 


IRDC308, 


.13- 


■1 


O.J \\J.\JZ.) 


Al ^ (C\ C\l\ 


Z, 


. 1 (_U.UU^ 


1. J 


d 1 

U. ID 


IRDC308, 


.13- 


■2 


Z. 1 (^U.Uj ) 


Al Q /^n C\l\ 


1 

1 , 


7 /^n di^ 


1 1 

1 . 1 


n 1 7 

U. i / 


IRDC308, 


.13- 


-3 


Z. 1 t^U.Uj ) 


Al s /^n n9~i 


1 

1 , 


S (d d9.\ 


1 n 
1 .u 


n 1 f\ 

U. iO 


IRDC309, 


,13- 


-1 


D.L 


A A A tr\ (M.^ 


a 

J, 


A (d Ad\ 


n 8 

u.o 


U.l J 


IRDC309. 


IS- 


-2 


Z.H- [\J.\H-} 


A< 9 {f\ C\f,\ 


J, 




u.o 


U. IH- 


IRDC309, 


IS- 


-3 


1 1 /A f\C\ 
1 . 1 (^U.Uj ) 


-'4-J.D (,U. 1 J ) 


■2 

J. 


A (d Ad\ 


U.J 


U.IO 


IRDC309, 


,37- 


-1 


IZ..Z. \ \J.l\J) 


Al 1 09'! 


z. 


'> (d 




n 41 

U.T-i 


IRDC309, 


.37- 


■2 


1 .U (^U.UJ ) 


-"4^.3 (^U.Uo ) 


z, 


1 (d 0(W 


U.J 


n 1 'X 

U. i J 


IRDC309, 


.37- 


■3 




-H-J.H \\J.\JZ,) 


z, 


1 (d df\\ 


1 c 

1 .J 


d 1 ^ 

U. i J 


IRDC310, 


.39- 


■1 


/ .y (U.Uj ) 


- jZ.O (,U.Ui ) 


J. 


A (d d'X\ 


9 1 

Z. 1 


n 1 ^ 

U. i J 


IRDC310, 


.39- 


■2 


j.u [yj.yjz.) 


- JZ. J (_U.Ui J 


z. 


(d d'W 


1 -O 


n 1 ^ 

U. i J 


IRDC312, 


.36- 


■1 


4. J (^U.UJ ) 


-JU.o (,U.Ui ) 


z, 


1 (d dA\ 


1 - J 


n 1 ^ 

U. i J 


IRDC312, 


.36- 


-2 


i.u \\J.\J^-} 


AQ A (Ci 


J, 




n 9 
u.z 


n 1 9 
u. iz 


IRDC313. 


.72- 


-1 


Al ((\ no\ 
. / yJ.yjA ) 




z, 


,\j \ yj,\j\j) 


9 9 
z.z 


n 19 
u.iz 


IRDC313, 


.72- 


-2 


J.J I^U.UJ ) 


AA A fC\ (^^\ 


9 

z, 


< en d'W 


9 n 

Z.U 


n 19 

U. IZ 


IRDC313 


.72- 


-3 


D.J \\J.\JA) 


A'i o ni \ 


z, 


.J \\J.\JZ.) 


9 

Z.J 


n 10 
u.iz 


IRDC313, 


.72- 


■4 


A 1 /^n 09^1 
^.1 t^U.UZj 


-'-tH.y.j \\.).\.} i ) 


z, 


d (d dd\ 


1 Q 

1 .y 


d 1 9 
U. 1 z 


IRDC316, 


.72- 


■1 


io.o l^U.U:> ) 


-jo.o [U.UL ) 


/I 


9 {d dl\ 
.Z \ \j.\.}Z.) 


J.D 


n 1 9 
u. iz 


IRDC316, 


.76- 


■1 


DA. 1 \\).\)D ) 


'jy.o \\.).\.) i ) 


A 




7 ^ 
/ .J 


n 1 9 

U. 1 z 


IRDC316, 


.72- 


■2 






a 

j>, 


f\ (d dX\ 


9 zl 
Z.^ 


n 1 9 

U. iZ 


IRDC316, 


.76- 


■2 


1 Q Zl iTl (\A \ 


Ad A (Ci ri] ^ 


/I 
^, 


1 /^n di\ 


A A 


n 1 ^ 

U. 1 J 


IRDC317, 


.71- 


-1 


H-.J (^U.UZ_; 


-'4Z. J (^U.Ul ^ 


z, 




9 n 

Z.U 


u.l J 


IRDC317. 


.71- 


-2 




-H-z. / t^u.uj. 


1 

i, 


■o i_u.uu_; 


j.i 


n 19 

U. IZ 


IRDC317. 


.71- 


-3 


\ A(C\ (Y)\ 


/to n r'O 0*^^ 


1 

1, 


■ O i^U.UU^ 


u. / 


n 1"^ 

u.l J 


IRDC320, 


.27- 


-1 


A (0 f)'X\ 

Z.H [V.Vj ) 


-JZ.i ^^U.UJ^ 


9 

z, 




u.o 


n 1"^ 

u.l J 


IRDC320, 


.27- 


-2 


\}.y (^U.UZ ) 


^9 ^ /^n CiA^ 


1 

1 , 


1 fd ^d'\ 


d A 


d 1 ^ 

U. 1 J 


IRDC320, 


.27- 


-3 




^9 (Ci 07^1 
-JjZ.O t^U.U / J 


1 

1 , 


^ /^n 9zL~i 
.J ^U.ZH-j 


d 1 


n 1 1 

VJ. L J 


IRDC321, 


.73- 


■1 


D.D t^U.UJ) ) 




z, 


^ /^n di\ 


9 1 
Z.J 


n 1 ^ 

U. i J 


IRDC321, 


.73- 


■2 






3 




1.8 


0.13 


IRDC321, 


.73- 


■3 


8.3 (0.03) 


-31.7 (0.01) 


3, 


.7 (0.03) 


2.1 


0.13 


IRDC013, 


.90- 


-1 


5.4 (0.04) 


23.1 (0.02) 


2, 


.4 (0.05) 


2.0 


0.12 


IRDC013, 


.90- 


-2 


1.8 (0.03) 


23.1 (0.05) 


3, 


.1(0.11) 


0.5 


0.12 


IRDC316. 


.45- 


-1 


4.5 (0.03) 


-43.5 (0.02) 


3, 


.3 (0.05) 


1.2 


0.15 


IRDC316, 


.45- 


-2 


4.7 (0.03) 


-43.3 (0.02) 


3, 


.2 (0.05) 


1.3 


0.14 


IRDC318, 


.15- 


-1 


3.7 (0.02) 


-42.7 (0.01) 


2, 


.5 (0.04) 


1.4 


0.13 


IRDC318, 


.15- 


-2 


3.3 (0.02) 


-42.0 (0.01) 


2, 


.1 (0.04) 


1.4 


0.13 


IRDC309, 


.94- 


-1 


8.5 (0.05) 


-59.5 (0.02) 


4, 


.2 (0.05) 


1.8 


0.16 



Notes: Columns are name, integrated area obtained from the gaussian fits, LSR velocity, full linewidth at half maximum, main beam brightness 
temperature, Icr rms values. 
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Table C.5. HCO"^ line parameters. 



Name 








AV 


T mb 


Icr rms 








(Kkm/s) 


(km/s) 




(km/s) 


(K) 


(K) 




1 2 


i 


A 1 (Ci CiA^ 




z 


7 


(0.07) 


1 c 

1 .J 





17 


TP r^r^'^ns 


1 ^ 


Z 


:>.^ (^U.Uj ) 


-'4 / .y \\)X)L ) 


z 


1 


(0.05) 


1 c 

1 .J 





16 




1 ^ 






zLS n (Ci n7~i 

-H-O.U t^U.UZj 


z 


3 


(0.07) 


1 7 

1 .z 





16 


TP r^r^aoQ 


1 'i 


1 
1 


4 wn 07^ 




Q 


1 


(0.41) 


n A 





16 




1 

i J 


7 

Z 


D,L yj.yjj ) 


-^D.Z [\J.\Jy} 


< 

J 


6 


(0.20) 


yj.D 





14 




1 

i J 






Al "X /"n 


Z 


6 


(0.14) 


n ^ 

V.D 





15 




J / 


1 
1 


1 ^ Q /Tl \A\ 


Al A 07 


z 


8 


(0.08) 


A ^ 





31 




J / 


Z 


1 .o (^U.Uj ) 


A'i fit /n 

-^■J.O yJXJD ) 


z 


8 


(0.13) 


u.o 





14 




J 1 


D 


J .y t^U.U J ) 


A'^ 7 (Ci 07^1 
-H-J. / (^U.UZ ) 


z 


2 


(0.06) 


7 A 





13 


IKUL^J iU 


3V 


1 
i 


O.J (^U.Uj ) 


n /^n 07^ 


z 


3 


(0.00) 







12 


Tpr^r'a 1 n 
iKUV_j iU 


Jy 


Z 


0.'+ (^U.UZJ 


-JZ.O (^U.Ui J 


z 





(0.00) 


^ n 

J.U 





12 


Tpp^/^a 1 n 


OD 


1 
1 


J. / (^U.Uj ) 


-J 1 .U (,U.Ui ) 


z 


7 


(0.04) 


1 Q 

1 .y 





12 


TDr^r"^ 1 


:)0 


Z 


1 1 /n fM \ 
i.i ^^U.LW-^ 


CI A (Ti 7 1 ^ 


c 
J 


7 


(0.44) 


n 1 





12 


Tpr^z-^a 1 a 


/Z 


1 
1 


7 7 /"n nA\ 
/.z t^u.uo^ 




z 


5 


(0.00) 


7 7 





11 




70 
/Z 


7 

Z 


7 1 /A (\A\ 




z 


5 


(0.00) 


7 7 
Z. / 





13 


TPr*r''5i 


70 

/Z 


:j 


Q 1 /n 

O.l \\J\}^) 


AA 7 fc\ ^^^\ 


z 


6 


(0.03) 


7 ft 
Z.o 





12 


TPT^r'a 1 a 


77 

/ Z 


A 


1 V\ (C\ C\X\ 
1 X) \\).\jD ) 




z 


5 


(0.00) 


7 ^^ 
z.o 





12 


Tpr\r'-2 1 fi. 

iKJjV_j iO 


77 
/ Z 


1 
i 




-JO. 1 \y)X)i ) 







(0.00) 


A f\ 





12 




7/^ 
/D 


1 




'jyX) \\.)X) i ) 


/I 





(0.00) 


7 n 





12 


Tpr^r'a 1 £. 
iKUV_j iO 


77 

/ Z 


Z 




-Jo.y \\.)X) L ) 


J 


4 


(0.05) 


7 7 
z.z 





12 




7/^ 


7 
z 


7^^ 7 iTl (\A \ 
Zj.Z l^U.UH-y 


-jy.j \\)X)L ) 


D 


5 


(0.03) 


J.J 





12 


Tpr^nQ 1 1 


71 
/ 1 


1 
1 


H-.J (^U.UZ_; 


-'41.7 (^U.Ul ) 


z 


1 


(0.00) 


7 n 

Z.u 





13 


Tpr\/^a 1 T 


71 


7 

Z 


J. < en ^7^ 


AO 1 (Tl 01^ 


z 


1 


(0.00) 


7 
Z.U 





12 


Tpr\/^ai 7 


71 


D 




47 7 ((\ ^7^ 


z 


1 


(0.00) 


1 n 

i.U 





12 




77 

Z / 


1 
1 


1 WO CWk 


-ji.U yjxj^) 


1 

1 


5 


(0.00) 


n 7 





12 


TPnr"^on 

IKJJL-jZU 


77 
Z / 


Z 


U. / (^U.uzj 


-DDXJ yJXJJ ) 


1 

1 


3 


(0.10) 


V.J. J 





12 


IKUV^ jZU 


77 

Z / 






^9 Q /A A9\ 


n 
u 


7 


(0.06) 


n s 





14 




73 


]^ 




^7 7 /^n nn 

- jZ. / ^U.Ui j 


z 


4 


(0.02) 


j.^^ 





13 


IRDC321 


73 


2 






5 





(0.06) 


1.4 





12 


IRDC321 


73 


3 


4.3 (0.10) 


-30.6 (0.05) 


1 


5 


(0.00) 


2.7 





14 


IRDC013 


90 


1 


5.4 (0.03) 


23.1 (0.01) 


2 


4 


(0.04) 


2.1 





45 


IRDC013 


90 


2 


2.1 (0.04) 


23.5 (0.07) 


3 


4 


(0.12) 


0.5 





20 


IRDC316 


45 


1 


4.7 (0.04) 


-43.2 (0.03) 


4 


3 


(0.09) 


1.0 





12 


IRDC316 


45 


2 


3.6 (0.03) 


-43.3 (0.04) 


4 


1 


(0.09) 


0.8 





13 


IRDC318 


15 


1 


1.9 (0.02) 


-43.3 (0.02) 


1 


8 


(0.00) 


0.9 





12 


IRDC318 


15 


2 


0.9 (0.02) 


-43.2 (0.06) 


1 


8 


(0.00) 


0.5 





12 


IRDC309 


94 


1 


9.9 (0.06) 


-60.3 (0.03) 


5 


4 


(0.08) 


1.7 





15 



Notes: Columns are name, integrated area obtained from the gaussian fits, LSR velocity, full linewidth at half maximum, main beam brightness 
temperature, Icr rms values. 



Table C.6. HNCO line parameters. 



Name 


/ Tmb dv 




A V 


T 


Icr rm 




(K km/s) 


(km/s) 


(km/s) 


(K) 


(K) 


IRDC309.37-1 


1.3 (0.07) 


-42.4 (0.12) 


2.3 (0.36) 


0.5 


0.15 


IRDC3 10.39-1 


0.3 (0.02) 


-53.8 (0.08) 


1.3 (0.20) 


0.2 


0.07 


IRDC312.36-1 


0.4 (0.02) 


-50.4(0.14) 


2.5 (0.26) 


0.1 


0.06 


IRDC313.72-1 


0.5 (0.02) 


-44.3 (0.06) 


1.4(0.19) 


0.3 


0.06 


IRDC313.72-2 


0.6 (0.02) 


-44.1 (0.13) 


2.9 (0.30) 


0.1 


0.05 


IRDC313.72-3 


0.3 (0.02) 


-43.6(0.11) 


2.0 (0.36) 


0.1 


0.06 


IRDC316.72-1 


0.9 (0.03) 


-39.1 (0.08) 


2.6(0.18) 


0.3 


0.06 


IRDC3 16.76-1 


1.0 (0.03) 


-39.7 (0.10) 


3.3 (0.29) 


0.2 


0.06 


IRDC317.71-1 


0.6 (0.03) 


-43.5 (0.09) 


2.0 (0.23) 


0.3 


0.07 


IRDC317.71-2 


1.1 (0.04) 


-43.7 (0.10) 


3.4 (0.31) 


0.3 


0.06 


IRDC320.27-1 


0.3 (0.01) 


-31.7 (0.07) 


1.3 (0.17) 


0.2 


0.06 


IRDC3 16.45-1 


0.3 (0.02) 


-43.4 (0.08) 


1.4 (0.19) 


0.2 


0.07 


IRDC3 16.45-2 


0.4 (0.02) 


-43.2 (0.07) 


1.5 (0.16) 


0.3 


0.07 


IRDC309.94-1 


1.1 (0.03) 


-59.4 (0.08) 


2.6 (0.22) 


0.3 


0.07 



Notes: Columns are name, integrated area obtained from the gaussian fits, LSR velocity, full linewidth at half maximum, main beam brightness 
temperature, Icr rms values. 
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Table C.7. SiO line parameters. 



Name 


ITuo dv 

1 IVl £} 




A V 




Icr rms 




(Kkm/s) 


(km/s) 


(km/s) 


(K) 


(K) 


IRDC309.37-1 


0.9 (0.08) 


-42.6 (0.2) 


2.8 (0.5) 


0.3 


0.18 


IRDC3 10.39-1 


0.6 (0.03) 


-53.3 (0.2) 


3.9 (0.4) 


0.1 


0.08 


1RDC3 12.36-1 


0.2 (0.03) 


-50.8 (0.2) 


2.0 (0.5) 


0.1 


0.08 


^ 79 1 

llvl_/V^ J ID, / Z- 1 




-^j.y [\J, i ) 


H-. J yJ.D ) 


yj.j 


U.Uo 


IRDC313.72-2 


1.7 (0.05) 


-AAA (0.1) 


5.1 (0.4) 


0.3 


0.08 


IRDC313.72-3 


1.3 (0.06) 


-41.9 (0.3) 


8.2 (0.8) 


0.1 


0.08 


IRDC313.72-4 


0.5 (0.03) 


-44.1 (0.2) 


3.4 (0.4) 


0.1 


0.08 


IRDC3 16.76-1 


2.4 (0.07) 


-39.5 (0.2) 


8.9 (0.7) 


0.2 


0.08 


1RDC317.71-3 


0.3 (0.03) 


-42.1 (0.2) 


1.9 (0.5) 


0.1 


0.08 


1RDC321.73-1 


2.1 (0.05) 


-32.1 (0.1) 


4.8 (0.3) 


0.4 


0.08 


IRDC309.94-1 


2.4 (0.05) 


-59.6 (0.2) 


5.0 (0.0) 


0.4 


0.10 



Notes: Columns are name, integrated area obtained from the gaussian fits, LSR velocity, full linewidth at half maximum, main beam brightness 
temperature, Icr rms values. 



Table C.8. H^CQ-" Une parameters. 



Name 


1 T\Aii dv 

J -* MB ' 


' Isr 




A V 


T 




f 


Icr rms 


6Y 






(km/s) 
















IRDC308.13-1 


0.6 (0.06) 


-46.7 (0.28) 


2, 


.5 (0.41) 


0.2 


0. 


,07 


0.08 


-0.8 (0.12) 


IRDC309.13-1 


0.8 (0.05) 


-44.7(0.14) 


2, 


.1 (0.31) 


0.3 


0. 


,20 


0.07 


-0.2 (0.09) 


IRDC309.13-3 


0.5 (0.05) 


-45.2 (0.12) 


1, 


.5 (0.39) 


0.3 


0. 


,17 


0.08 


-2.1 (0.35) 


Tr>T~\/~''3no 11 1 
IKUCJUy.J /-I 


z.o (U. 13) 


-42.0 (V.W) 


z, 


. 1 (U./O) 


1 1 
i.i 


A 

u. 


,zu 


U.lD 


A Q /A A/:\ 

-U.o (U.Uo) 


IRDC309.37-3 


1.3 (0.06) 


-42.4 (0.15) 


3, 


.2 (0.34) 


0.3 


0. 


,10 


0.08 


-1.2 (0.06) 


IRDC310.39-1 


1.9 (0.04) 


-52.4 (0.04) 


1 


.9 (0.11) 


0.9 


0. 


19 


0.08 


-0.6 (0.03) 


1RDC3 10.39-2 


1.2 (0.05) 


-52.3 (0.09) 


2, 


.2 (0.20) 


0.5 


0. 


12 


0.08 


-0.4 (0.04) 


1RDC312.36-1 


1.3 (0.06) 


-50.6(0.13) 


3, 


.0 (0.35) 


0.4 


0. 


12 


0.08 


-0.3 (0.04) 


1RDC3 12.36-2 


0.2 (0.03) 


-48.8 (0.20) 


1, 


.4 (0.38) 


0.1 


0. 


11 


0.08 


-2.5 (0.49) 


IRDC3 13.72-1 


1.1 (0.08) 


-43.9 (0.19) 


2, 


.5 (0.63) 


0.4 


0. 


10 


0.08 


-0.8 (0.11) 


IRDC313.72-2 


1.1 (0.05) 


-43.9(0.10) 


2, 


.2 (0.31) 


0.4 


0. 


,11 


0.08 


-0.8 (0.06) 


IRDC313.72-3 


0.7 (0.04) 


-43.5(0.10) 


1, 


.8 (0.21) 


0.4 


0. 


,09 


0.08 


-0.7 (0.07) 


IRDC313.72-4 


0.8 (0.05) 


-43.3(0.12) 


2, 


.2 (0.30) 


0.3 


0. 


,09 


0.08 


-1.2 (0.09) 


IRDC3 16.72-1 


4.4 (0.06) 


-38.8 (0.04) 


3, 


.1 (0.10) 


1.2 


0. 


,22 


0.08 


0.1 (0.01) 


IRDC3 16.76-1 


5.2 (0.06) 


-39.7 (0.03) 


3, 


.2 (0.09) 


1.5 


0. 


18 


0.08 


0.7 (0.01) 


IRDC3 16.72-2 


2.7 (0.05) 


-39.5 (0.04) 


2, 


.2 (0.09) 


1.1 


0. 


,33 


0.08 


0.6 (0.02) 


1RDC3 16.76-2 


3.1 (0.05) 


-40.5 (0.04) 


2, 


.4 (0.10) 


1.1 


0, 


,26 


0.08 


1.1 (0.02) 


1RDC317.71-1 


1.3 (0.04) 


-43.3 (0.07) 


1, 


.9 (0.17) 


0.6 


0. 


18 


0.08 


1.3 (0.07) 


1RDC317.71-2 


1.7 (0.05) 


-43.8 (0.05) 


1, 


.9 (0.16) 


0.8 


0. 


,23 


0.08 


1.7 (0.07) 


IRDC317.71-3 


0.4 (0.04) 


-43.7 (0.20) 


1, 


.7 (0.47) 


0.2 


0. 


,08 


0.09 


1.0 (0.19) 


IRDC320.27-1 


0.7 (0.03) 


-31.9 (0.05) 


1, 


.1 (0.12) 


0.5 


0. 


,28 


0.09 


0.8 (0.10) 


IRDC321.73-1 


1.6 (0.04) 


-32.1 (0.04) 


1, 


.5 (0.09) 


0.9 


0. 


,20 


0.08 


-0.5 (0.03) 


IRDC321.73-2 


0.6 (0.04) 


-33.4 (0.08) 


1, 


.3 (0.22) 


0.4 


0. 


,16 


0.08 


0.7 (0.10) 


IRDC321.73-3 


1.2 (0.06) 


-31.6(0.14) 


2, 


.5 (0.31) 


0.4 


0. 


10 


0.09 


1.0 (0.07) 


1RDC013.90-1 


1.0 (0.04) 


23.3 (0.09) 


1, 


.9 (0.18) 


0.5 


0. 


14 


0.08 


-0.2 (0.04) 


1RDC013.90-2 


0.3 (0.02) 


22.3 (0.03) 


0, 


.4 (0.10) 


0.6 


0. 


,35 


0.08 


1.2 (0.68) 


1RDC3 16.45-1 


0.8 (0.05) 


-43.5 (0.12) 


1, 


.9 (0.25) 


0.4 


0. 


16 


0.08 


0.2 (0.06) 


IRDC3 16.45-2 


1.1 (0.05) 


-43.3 (0.10) 


2, 


.2 (0.23) 


0.4 


0. 


,21 


0.08 


-0.0 (0.05) 


IRDC318.15-1 


0.7 (0.04) 


-42.5 (0.09) 


1, 


.6 (0.21) 


0.4 


0. 


,17 


0.08 


-0.7 (0.08) 


IRDC318.15-2 


0.6 (0.04) 


-42.0(0.10) 


1, 


.5 (0.20) 


0.3 


0. 


,19 


0.08 


-1.1 (0.13) 


IRDC309.94-1 


1.7 (0.06) 


-59.3 (0.07) 


2, 


.1 (0.17) 


0.7 


0. 


,25 


0.11 


-0.9 (0.05) 



Notes: Columns are name, integrated area obtained from the gaussian fits, LSR velocity, full linewidth at half maximum, main beam brightness 
temperature, optical depth, Icr rms values, SV = (V//co+ - V„i3co+)/^Vh"co->-- 
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Table C.9. Integrated areas, measured by summing the channels between suitable velocity limits under the line. 



Name 


NiH+ 




HC,N 


HNC 


HCO+ 


HCN 


HNCO 




SiO 


H'^CN 


CH^&H 




( K km/s) 


(K kni/s) 


(K km/s) 


(K kni/s) 


(K km/s) 


(K km/s) 


(K km/s) 


(K kni/s) 


(K km/s) 


(K km/s) 


(K km/s) 


IRDC308.13-1 


4.2 






4.8 


6.2 


5.3 












IRDC308.13-2 


2.1 






2.6 


4.2 


2.9 




0.6 








IRDC308.13-3 


1.6 






2.9 


3.7 


3.3 












IRDC309.13-1 


5.0 






3.1 


4.1 


4.5 




0.4 








IRDC309.13-2 


2.8 






2.5 


3.2 


3.9 












IRDC309.13-3 


1.5 






0.8 


1.9 


2.3 












IRDC309.37-1 


29.1 




1.6 


13.6 


17.9 


15.4 


1.4 


3.2 








IRDC309.37-2 








1.2 


2.4 


2.4 












IRDC309.37-3 


4.6 






4.8 


6.8 


6.0 












IRDC3 10.39-1 


13.4 




1.9 


8.1 


13.3 


12.5 


0.2 


1.7 


0.4 






IRDC3 10.39-2 


5.4 




1.2 


5.5 


8.7 


9.0 




1.4 








IRDC3 12.36-1 


6.2 




0.2 


4.9 


6.6 


6.6 


0.6 


0.7 








IRDC3 12.36-2 


1.3 






1.9 


1.0 


2.1 












IRDC3 13.72-1 


8.0 




1.7 


6.9 


11.3 


7.8 


0.6 


1.1 


1.6 






IRDC313.72-2 


8.6 


0.3 


1.2 


6.7 


10.6 


8.0 


0.5 


0.9 


1.9 






IRDC313.72-3 


7.1 




0.6 


7.3 


9.1 


7.3 


0.4 


0.6 


1.3 






IRDC313.72-4 


5.9 






5.7 


9.3 


5.9 






0.6 






IRDC3 16.72-1 


18.4 




1.1 


17.1 


23.2 


18.6 


0.8 


2.8 




1.9 


0.3 


IRDC3 16.76-1 


50.7 


0.9 


4.9 


34.1 


34.6 


33.8 


1.1 


7.1 


2.3 


4.8 


0.7 


IRDC3 16.72-2 


11.8 




0.8 


10.2 


14.1 


10.9 




1.9 






0.3 


IRDC3 16.76-2 


27.4 


0.5 


2.9 


21.4 


23.3 


27.4 




7.1 




3.3 


0.3 


IRDC317.71-1 


10.5 




0.4 


5.0 


5.6 


2.7 


0.7 


0.8 








IRDC3 17.7 1-2 


21.9 




2.9 


8.1 


5.6 


4.9 


1.3 


1.5 




1.5 


0.57 


IRDC317.71-3 


3.5 






2.3 


2.5 


1.8 






1.0 






IRDC320.27-1 


2.5 






2.6 


2.6 


3.8 












IRDC320.27-2 


1.7 






1.1 


0.7 


1.7 












IRDC320.27-3 


0.6 






0.8 


0.3 


0.8 












IRDC321.73-1 


11.0 




1.3 


7.0 


10.5 


11.1 




1.4 


2.5 






IRDC321.73-2 


5.5 




1.7 


6.9 


8.4 


7.9 




0.9 








IRDC321.73-3 


7.5 




1.3 


8.4 


10.4 


8.9 




0.8 








IRDC013.90-1 


5.6 




0.7 


5.6 


6.3 


5.4 




0.9 








IRDC013.90-2 


2.4 






1.4 


0.8 


1.9 












IRDC316.45-1 


7.0 




1.0 


4.7 


5.1 


4.0 


0.1 


0.7 








IRDC3 16.45-2 


6.2 






4.9 


3.4 


4.0 


0.4 


0.7 








IRDC318.15-1 


3.7 






3.9 


3.0 


3.3 




0.7 








IRDC318.15-2 


4.6 






3.9 


1.8 


2.9 




0.5 








IRDC309.94-1 


22.0 




3.7 


9.4 


11.0 


11.8 


0.9 


2.4 




2.1 





Notes: Columns are name, integrated area for every species. 



